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We report measurements of the exclusive cross section for e + e~ —* D°D*~tt + as a function 
of center-of-mass energy from the D°D*~ir + threshold to 5.2 GeV with initial-state radiation. No 
evidence is found for F(4260) — > D°D*~tt + decays. The analysis is based on a data sample collected 
with the Belle detector at or near a center-of-mass energy of 10.58 GeV with an integrated luminosity 
of 695 fb" 1 at the KEKB asymmetric-energy e + e~ collider. 
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Studies of exclusive open charm production near 
threshold in e + e~ annihilation provide important infor- 
mation on the dynamics of charm quarks and on the 
properties of the ^ states. During the past three years 
numerous measurements of exclusive e + e~ cross sections 
for charmed hadron pairs have been reported. Most of 
these measurements were performed at i?-factories using 
initial-state radiation (ISR). Belle presented the first re- 
sults on the e + e~ cross sections to the DD, D + D*~ 
D* + D*~ , D Q D~ir + (including the first observation of 
tp(U15Y -> DD* (2460) decays) and A+A~ final 

states Q- BaBar measured e + e~ cross sections to DD 
and recently to the DD*, D*D* final states @,0]. CLEO- 
c performed a scan over the energy range from 3.97 to 
4.26 GeV and measured exclusive cross sections for the 
DD, DD* and D*D* final states at thirteen points with 
high accuracy • The measured open charm final states 
nearly saturate the total cross section for charm hadron 
production in e + e~ annihilation in the ^fs region up to 
~ 4.3 GeV. In the energy range above ~ 4.3 GeV some 
room for contributions to the -0(4415) state from unmea- 
sured channels still remains. The exclusive cross sections 
for charm strange meson pairs have been measured to be 
an order of magnitude smaller than charm meson pro- 
duction Q. Charm baryon-antibaryon pair production 
occurs at energies above 4.5 GeV. 

Another motivation for studying exclusive open charm 
production is the existence of a mysterious family of 
charmonium-like states with masses above open-charm 
threshold and quantum numbers J PC = 1 . Although 
these have been known for over four years, the nature 
of these states, found in e + e~ — > 7r + 7r~ J/ip(ip(2S))^isR 
processes, remains unclear. Among them are the Y (4260) 



F(4008) and the F(4660) seen by Belle [HEl. 

No clear evidence for open charm production associ- 
ated with any of these states has been observed. In fact 
the F(4260) peak position appears to be close to a local 
minimum of both the total hadronic cross section [16JI and 
of the exclusive cross section for e + e~ — > D*D* 0, 0]. 
The AT(4630), recently found in the e + e~ — ► A+Aj cross 
section as a near-threshold enhancement Q, has a mass 
and width (assuming the AT (4630) to be a resonance) 
consistent within errors with those of the F(4660), sup- 
porting explanation that the XU630) is F(4660) 
or 0(25)/ o (98O) bound state However, this co- 

incidence does not exclude other interpretations of the 
X (4630), for example, as a conventional charmonium 
state [l!| or as baryon-antibaryon threshold effect [2fl |. 
point-like baryons [2l[ , or as a tetraquark state [22j] • 

The absence of open charm decay channels for Y 
states, large partial widths for decay channels to char- 
monium plus light hadrons and the lack of available 
J PC = 1 charmonium levels are inconsistent with the 
interpretation of the Y states as conventional charmonia. 
To explain the observed peaks, some models assign the 
3 3 L»i(4350), 5 3 Si(4660) with shifted masses [H, other 
explore c oup led-channel effects and rescattering of charm 
mesons US- More exotic suggestions include hadro- 



charmonium |25lk multiquark states, such as a 
tetraquark [26| and DD\ or D°D*° molecules 



27] 



state observed by BaBar 0,113, confirmed by CLEO 
H3 and Belle the y(4350) discovered by BaBar 
and confirmed by Belle (l5j : and two structures, the 



11, 

14] 



cq] [cq] 
One 

of the most popular exotic options for the Y states are 
the hybrids expected by LQCD in the mass range from 
4.2 — 5.0GeV/c 2 [28j]. In this context, some authors ex- 
pect the dominant decay channels of the Y(4260) to be 
^(4260) -► DW^Wtt. 

In this paper we report a measurement of the exclu- 
sive e + e~ — > D°D*~tt + cross section as a function of 
center-of-mass energy from the D°D*~tt + threshold to 
5.2 GeV, as part of our studies of the exclusive open- 
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charm production in this mass range. The analysis is 
based on a data sample collected with the Belle detec- 
tor [2§| at the T(4S f ) resonance and nearby continuum 
with an integrated luminosity of 695 fb -1 at the KEKB 
asymmetric-energy e + e~ collider [30| . 

We employ the reconstruction method that was used 
for e+e~ — > DD and e + e~ — ► D°D~ir + exclusive 
cross section measurements @, 0|- We select e+e~ — > 
D°D* _ 7r + 7isR signal events by reconstructing the D°, 
D*~ and n + mesons. In general the 7isr is not required 
to be detected: instead, its presence in the event is in- 
ferred from a peak at zero in the spectrum of recoil mass 
squared against the D°D*~tt + system. The square of the 
recoil mass is defined as: 

Ml con (D a D*-n+) = (E c . m . - E D o D ,- w+ f-pl OD ,. w+ (l) 

Here -E c .m. is the initial e + e _ center-of-mass (cm.) en- 
ergy, i?£)0£ ) *- 7r + and Pd°d*-tt+ are the cm. energy 
and momentum of the D°D*~ir + combination, respec- 
tively. To suppress backgrounds two cases are consid- 
ered: (1) the 7isr is outside of the detector acceptance 
and the polar angle for the D°D*~ir + combination in 
the cm. frame is required to be \cos(9]jo D *- 7T +)\ > 
0.9; (2) the fast 7isr is within the detector acceptance 
(|cos(0 D o£ ) »- 7r +)| < 0.9), in which case it is required to be 
detected and the mass of the D°D*~ ir + -)isn combination 
must be greater than (-E^m. — 0.58GcV/c 2 ). To sup- 
press background from e + e~ — > D°D*~ir + (2n)ir 7isr 
(n > 1) processes we exclude events that contain addi- 
tional charged tracks that are not used in D°, D*~ or ir + 
reconstruction. 

All charged tracks are required to originate from the 
vicinity of the interaction point (IP); we impose the re- 
quirements dr < 1 cm and |c?z| < 4cm, where dr and \dz\ 
are the impact parameters perpendicular to and along 
the beam direction with respect to the IP. Charged kaons 
are required to have a ratio of particle identification likeli- 
hood, Vk = £k/{£k+£tt), larger than 0.6 [3l|. Charged 
tracks not identified as kaons are assumed to be pions. 

Kg candidates are reconstructed from ir + ir~ pairs with 
an invariant mass within 10 MeV/c 2 of the Kg mass. The 
distance between the two pion tracks at the K s vertex 
must be less than 1 cm, the transverse flight distance from 
the IP is required to be greater than 0.1cm, and the 
angle between the Kg momentum direction and the flight 
direction in the x—y plane should be smaller than 0.1 rad. 

Photons are reconstructed from showers in the electro- 
magnetic calorimeter with energies greater than 50 McV 
that are not associated with charged tracks. ISR pho- 
ton candidates are required to have energies greater than 
2.5 GeV. Pairs of photons are combined to form 7r° can- 
didates. If the mass of a 77 pair lies within 15 MeV/c 2 of 
the 7T° mass, the pair is fitted with a ir° mass constraint 
and considered as a tt° candidate. 

D a candidates are reconstructed using five decay 
modes: K~ir + , K~K + , K~ir~ir + ir + , K s it + it~ and 
K~ir + ir°. A ±15McV/c 2 mass window is used for all 
modes except for K~~ir~ir + ir + where a ±10MeV/c 2 re- 



quirement is applied (~ 2.5 er in each case). D + candi- 
dates are reconstructed using K~n + ir + and K s n + decay 
modes [HI; a ±15MeV/c 2 mass window is used for both 
D + modes. To improve the momentum resolution of D 
meson candidates, final tracks are fitted to a common 
vertex with a mass constraint on the D° or D + mass. 
D* candidates are selected via the D* + — > D Q ir + and 
D*° — > D°tt° (for background study) decay modes with 
a ±2MeV/c 2 D* — D mass-difference window (~ 3 c). 
A mass- and vertex-constrained fit is also applied to D* 
candidates. 

To remove contributions from the e + e~ — > 
£)*+ £)*~7 ISR process, we exclude D°ir + combina- 
tions with invariant mass within ±5MeV/c 2 of the 
nominal D* + mass. 

D°, D + , D*~ and D* a mass sidebands are selected 
for the background study; these are four times as large 
as the signal region and are subdivided into windows 
of the same width as the signal. To avoid signal over- 
subtraction, the selected D sidebands are shifted by 
30McV/c 2 (20MeV/c 2 for the D° -> K--k--k+ti+ mode) 
from the signal region. The D candidates from these 
sidebands are refitted to the central mass value of each 
window. D* sidebands arc shifted by 4McV/c 2 to the 
higher mass side of the signal region. 

The distribution of Af 2 ocoil ( J D°D*-7r+) for the sig- 
nal region in the data for M D a D ,-^+ < 5.2GeV/c 2 
is shown in Fig. [T]a). A clear peak corresponding to 
the e + e~ — ► D°D*~n + -fisR process is evident around 
zero. The shoulder at positive values is due to e + e~ — > 
D D* _ 7t + 7isr + (n)7r° + (771)7 events. We define the 
signal region for Ai r 2 ccoil (Z3 £'* _ 7r + ) by a tight require- 
ment ±0.7( GeV/c 2 ) 2 around zero to suppress the tail 
from such events. The invariant-mass distribution of 
~tt + 7isr combinations in the data after the re- 
quirement on M 2 ecoil (D°_D*~7r + ) and the polar angle dis- 
tribution of D°D*~7T + combinations shown in Fig. [T]b), 
c) are typical of ISR production and are in agreement 
with the MC simulation. The M d q d *- v + spectrum ob- 
tained after all the requirements is shown in Fig. [2j 

The contribution of multiple entries after all the re- 
quirements is found to be less than 6%. In such case the 
single D°D*~tt + combination with the minimum value of 
Xtot = xL(D°) + Xm(d—) is chosen, where x\ [[D0) and 
Xm(d--) correspond to the mass fits for D° and D*~ 
candidates. 

The following sources of background are considered: 

(1) combinatorial background under the D°(D*~) 
peak combined with a real D*~(D°) coming from 
the signal or other processes; 

(2) both D° and D*~ are combinatorial; 

(3) the reflection from the processes e + e~ — * 
D°D*~7r + 7rJJ liss 7i S R, where the 7^^, is not recon- 
structed, including D*° — > D°7T^ liss decays; 

(4) the reflection from the process e + e~ — > 
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FIG. 1: a) The distribution of M? ecoil {D° D*~ -k + ). b) The 
mass spectrum of D°D* ~7t + 7isr combinations, c) The polar 
angle distribution of D°D*~tv + combinations. Histograms 
show the normalized M D o and M D ,~ sideband contributions. 
The selected signal windows are illustrated by vertical dotted 
lines. 
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FIG. 2: The M D o D ,- n + spectrum. The histogram shows the 
normalized M n o and M n *- sideband contributions. 



D* D*-7t+7isr, followed by D*° -> D°7, where 
the low-momentum 7 is not reconstructed; 



(5) ^ 



D°D* 7r + 7r where the energetic ir° is 



misidentified as a single 7isr- 

The contribution of background (1) is extracted using 
the D° and D*~ sidebands. Background (2) is present 
in both the M D o and M D *- sidebands and is, thus, sub- 
tracted twice. To take into account this over-subtraction 



we use a two-dimensional sideband region, when events 
are selected from both the M D o and the M D »- side- 
bands. The total contribution from the combinatorial 
backgrounds (1-2) is shown in Figs, [TJ [5] as a hatched 
histogram. 

Most of the background (3-4) events are suppressed 
by the tight requirement on Af 2 ocoil (D°Z}*~7r + ). The re- 
mainder of background (3) is estimated directly from the 
data by applying a similar reconstruction method to the 
isospin-conjugate process e + e~ — > -D + -D*~7r + 7r miss 7isR- 
Since there is a charge imbalance in the D + D*~tt + final 
state, only events with a missing extra 7r~ iss can con- 
tribute to the M 1 2 ocoil (£> + L>*~7r + ) signal window. To ex- 
tract the level of background (3), the D + D*~tt + mass 
spectrum is rescaled according to the ratio of D~ and 
D° reconstruction efficiencies and an isospin factor of 
1/2. A negligibly small contribution of background (3) 
is found: only one event with M D + D ,- n + < 5.2GeV/c 2 . 
Uncertainties in this estimate are included in the system- 
atic error. The remainder of background (4) is estimated 
from the data assuming isospin symmetry. We measure 

■4 D%°) by 
Only three 



V!.() 



the process e + e~ — > L>* + I?* 7r + 7 ISR (D 
applying a similar reconstruction method, 
events with M D „ + yj„o 7r+ < 5.2GeV/c 2 are found in the 
data. Thus the contribution of background (4) is also 
found to be negligibly small; uncertainties in this esti- 
mate are included in the systematic error. 

The contribution of background (5) is determined from 
the data using fully reconstructed e + e~ — > D a D*~ir + Tr° 
events including the reconstruction of an energetic 7r°. 
Only one event with M D o D *- 7T + < 5.2GeV/c 2 and 
Mdod"-h+-k° > 10GcV/c 2 is found in the data. As- 
suming a uniform ir° polar angle distribution, this back- 
ground contribution to the \cos(0 D o D „- 7T +)\ > 0.9 signal 
sub-sample (case 1) is 1 event/Q^o ~ 0.2 events in the 
entire Mdo£>*- w + mass range, where rj^o is the ir° recon- 
struction efficiency. The probability of 7r° — > 7 misiden- 
tification due to asymmetric 7r° — > 77 decays is also esti- 
mated to be C 1. Thus the contribution of background 
(5) is found to be negligibly small; uncertainties in this 
estimate are included in the systematic error. 

The e + e~ — > D°D*~tt + cross section is extracted from 
the background subtracted D°D*~tt + mass distribution 

<j(e+e- -► D°D* 



dN/dm 
r] tot dL/dm 



(2) 



where m = Mjyout- 7T + , dN/dm is the mass spectrum 
obtained without corrections for resolution and higher- 
order radiation, 7j to t is the total efficiency, and the factor 
dL/dm is the differential ISR luminosity 33]. The total 
efficiency determined by MC simulation grows quadrati- 
cally with energy from 0.007% near threshold to 0.036% 
at 5.2 GeV/c 2 . The resulting e+e" -> D°D*-tt+ exclu- 
sive cross section averaged over the bin width is shown 
in Fig. [3] with statistical uncertainties only. Since the 
bin width is much larger than the M D a D ,-^+ resolu- 
tion, which varies from ~ 3McV/c 2 around threshold to 



6McV/c 2 at M D o D . 



5.2GcV/c , no correction 
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for resolution is applied. 

The systematic errors for the cr(e + e~ — > D° D*~ir + ) 
measurement are summarized in Table HI The system- 

TABLE I: Contributions to the systematic error on the 
e + e~ — > D°D*~ir + cross section. 



Source 


Error, [%] 


Background subtraction 


±3 


Cross section calculation 


±6 


B(D) 


±3 


Reconstruction 


±7 


Kaon identification 


±2 


Total 


±10 



atic errors associated with the background (1-2) subtrac- 
tion are estimated to be 2% due to the uncertainty in the 
scaling factors for the sideband subtractions. This is es- 
timated from fits to the Mpo and M D *- distributions 
in the data that use different signal and background pa- 
rameterizations. Uncertainties in backgrounds (3-5) are 
estimated conservatively to be each smaller than 1% of 
the signal. The systematic error ascribed to the cross 
section calculation includes a 1.5% error on the differ- 
ential luminosity and a 6% error in the total efficiency 
function. Another source of systematic errors is the un- 
certainties in track and photon reconstruction efficiencies 
(1% per track, 1.5% per photon and 5% per K s ). Other 
contributions come from the uncertainty in the identifi- 
cation efficiency and the absolute D° and D*~ branching 
fractions [34[. The total systematic uncertainty is 10%. 

We perform a likelihood fit to the M D o D *-^+ distri- 
bution where we parameterize a possible -0(4415) sig- 
nal contribution by an s-wave relativistic Breit-Wigncr 
(RBW) function with a free normalization. We use PDG 
values [Hj to fix its mass and total width. To take a 
non-resonant D°D*~tt + contribution into account we use 
a threshold function ^ ' AI — m D o — m D *~ — with a 
free normalization. Finally, the sum of the signal and 
non-resonant functions is multiplied by a mass-dependent 
second-order polynomial efficiency function and differen- 
tial ISR luminosity. 

The fit yields 14.4 ± 6.2(stat.)^g'°(sys.) signal events 
for the 0(4415) state. The statistical significance for the 
■0(4415) signal is determined to be 3.1a from the quantity 
— 2 ln(£o/£max), where £ ma x is the maximum likelihood 
returned by the fit, and Cq is the likelihood with the am- 
plitude of the RBW function set to zero. The goodness 
of the fit is \ 2 /n.d.f = 1.17. The systematic errors of 
the fit yield are obtained by varying the mass and to- 
tal width within their uncertainties, histogram bin size 
and the parameterization of the background function and 
efficiency. 

We calculate the peak cross section for the e + e~ — > 
-0(4415) — > D°D*~n + process at E c , m , = to^( 4415 ) from 
the amplitude of the RBW function in the fit to be 



cr(e+e- -> 0(4415)) x £(0(4415) -> D°D*-tt+) < 0.76 
nb at the 90% C.L. Using a{e + e" -> 0(4415)) = 
127r/77i^ (4415) x (B cc ) and PDG values of the 0(4415) 
mass, full width and electron width [34[ we found B ee x 
£(0(4415) -► D°D*-tt+) < 0.99 x 10~ 6 at the 90% C.L 
and £(004415) -► D°D*-tt + ) < 10.6% at the 90% C.L. 
All presented upper limit values include systematic un- 
certainties. For illustration we include the corresponding 
fit function on the cross section distribution plot shown 
in Fig. [3 
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FIG. 3: The exclusive cross section for e e~ — > D D*~n 
averaged over the bin width with statistical uncertainties only. 
The fit function corresponds to the upper limit on 0(4415) 
taking into account systematic uncertainties. The solid line 
represents the sum of the signal and threshold contributions. 
The threshold function is shown by the dashed line. 

To obtain limits on the decays X — > D°D*~ tt + , 
where X denotes F(4260), F(4350), F(4660) or AT (4630) 
states, we perform four likelihood fits to the M D a D *- jr + 
spectrum each with one of the X states, the -0(4415) 
state and a non-resonant contribution. For fit functions 
we use the sum of two s-wave relativistic RBW func- 
tions with a free normalization and a threshold func- 
tion \/M — mjjo — Tfiu*— — 171^+ with a free normaliza- 
tion. The sum of the signal and non-resonant functions 
is multiplied by the mass-dependent second-order poly- 
nomial efficiency function and differential ISR luminosity. 
For masses and total widths of the V(4260) and 0(4415) 
states we use PDG values [Hj]. The corresponding pa- 
rameters of the F(4660), y(4350) and X(4630) states 
are fixed from Ref . [f| [35| , respectively. 

The significances for the Y*(4260), Y(4350), F(4660) 
and AT(4630) signal are found to be 0.9er, 1.4er, O.ler 
and 1.8cr, respectively. The calculated upper limits (at 
the 90% C.L.) on the peak cross sections for e + e~ — > 
X — > D°D*~tt + processes at E c , m , = mx are presented 
in Table [TTJ Using fixed values of X masses and full 
widths we obtain upper limits on the B ee x B(X — > 
D°D*-ir+) at the 90% C.L. Finally, for the F(4260) 
state we estimate the upper limit on £>(V(4260) — > 
D°D*-tt+)/B(Y(4260) -> Tr+Tr-J/0) at the 90% C.L. 
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using B cc x r(7r+7r-J/V>) H. For the F(4350) and 
F(4660) states we calculate B(X -► £> £>*-7r+)/£(X -> 
Tt + Tr~tp(2S)) at the 90% C.L. taking into account B ee x 
r(n + n~tp(2S)) [H. All upper limits presented in Ta- 
ble [II] are determined by choosing the maximum signal 
amplitudes that emerge from: varying masses and widths 
of the X states within their uncertainties; varying the his- 
togram bin size; and changing the parameterizations of 
the background & efficiency functions. 

To estimate the effects of possible interference be- 
tween final states we also performed a fit to the 
^d°d*-tt+ spectrum that includes complete interference 
between the -0(4415) RBW amplitude and a non-resonant 
D°D*~ir + contribution. We found two solutions both 
with x 2 /n.d.f = 1.28; the interference is constructive for 
one solution and destructive for the other. From the fit 
with destructive interference we find an upper limit on 
the peak cross section for e + e~ — * 0(4415) — * D°D*~tt + 
process to be cr(e+e~ — > 0(4415)) x 6(0(4415) — * 
D D*-n+) < 1.93 nb at the 90% C.L. 

In addition we performed four likelihood fits to the 
M D o D ,~n+ spectrum with complete interference be- 
tween the X and -0(4415) states' RBW amplitudes 
and a non-resonant D°D*~ir + contribution. We found 
four solutions for each fit with similar goodncss-of-fit 
(X 2 / n.d.f =1.39, 1.23, 1.39 & 1.21) and obtained the 
upper limits on the peak cross sections for e + e~ — ► 
X -> D°D*-tt+ process to be a(e+e~ -> X) X B(X -> 



D°D*-tt + ) less than 1.44, 1.92, 1.38 and 0.98 nb at the 
90% C.L. for F(4260), F(4350), F(4660) and A(4630), 
respectively. 

In summary, we report the first measurement of the 
e + e~ — > D q D*~it + exclusive cross section over the 
center-of-mass energy range from 4.0 GcV to 5.2 GcV. 
We calculate an upper limit on the peak cross sec- 
tion for the e + e~ — ► 0(4415) — ► D°D*~tt + process 
at -Ecm. = ^v( 441 5) to k e 0.76 nb at the 90% C.L. 
The values of the amplitude of the F(4260), y(4350), 
y(4660) and A(4630) signal function obtained in the 
fit to the M D o D »-^+ spectrum are found to be con- 
sistent with zero within errors. We see no evidence 
for y(4260) — > D°D*~tt + decays as predicted by hy- 
brid models and obtain the upper limit S(F(4260) — > 
D D*-n + )/B(Y {4260) 7T+7T-J/0) < 9 at the 90% 
C.L. 
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